Spontaneous symmetry breaking 
in amnestically induced persistence 
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We investigate a recently proposed non-Markovian random walk model characterized by loss 
of memories of the recent past and amnestically induced persistence. We report numerical and 
analytical results showing the complete phase diagram, consisting of 4 phases, for this system: 
(i) classical nonpersistence, (ii) classical persistence (iii) log-periodic nonpersistence and (iv) log- 
periodic persistence driven by negative feedback. The first two phases possess continuous scale 
invariance symmetry, however log-periodicity breaks this symmetry. Instead, log-periodic motion 
satisfies discrete scale invariance symmetry, with complex rather than real fractal dimensions. We 
find for log-periodic persistence evidence not only of statistical but also of geometric self-similarity. 



Nonpersistent random walkers with negative feedback 
tend not to repeat past behavior ll| , but what happens 
when they forget their recent past [2|? Remarkably, they 
become persistent for sufficiently large memory losses. 
This recently reported phenomenon of amnestically in- 
duced persistence 0, Q allows log-periodic [1| superdif- 
fusion 0, 0, 0] driven by negative feedback. Its practical 
importance stems from the conceptual advance of quanti- 
tatively relating, on a causal level, two otherwise appar- 
ently unconnected phenomena: repetitive or persistent 
behavior on the one hand, and recent memory loss on 
the other, e.g., in Alzheimer's disease 0j. Precisely how 
does persistence depend on recent memory loss? Here, we 
answer this question and report numerical and analytical 
results showing the complete phase diagram for the prob- 
lem, comprising 4 phases: (i) classical nonpersistence (ii) 
classical persistence, (iii) log-periodic nonpersistence and 
(iv) log-periodic persistence driven by negative feedback. 
The former two phases possess continuous scale invari- 
ance symmetry, which breaks down in the other two. 

Random walkers without memory have a mean square 
displacement ( x 2 ) that scales with time t according to 
( x 2 ) ~ t 2H , with Hurst exponent H = 1/2 as demanded 
by the Central Limit Theorem, assuming finite moments. 
Hurst exponents H > 1/2 indicate persistence and can 
arise due to long-range memory. Most random walks with 
and without memory display continuous scale invariance 
symmetry, i.e., continuous scale transformations by a 
"zoom" factor A leave the Hurst exponent unchanged: 

t 2H ^ X 2H t 2H ag t Xt 

Schiitz and Trimper [l[ pioneered a novel approach for 
studying walks with long-range memory [fj, [7], |8|, [9( , which 
we have adapted [sj for studying memory loss. Consider 
a random walker that starts at the origin at time to = 0, 
with memory of the initial ft time steps of its complete 
history (0 < / < 1). At each time step the random 
walker moves either one step to the right or left. Let 



Vt = ±1 represent the "velocity" at time t, such that the 
position follows 



x t+1 = x t + v t +i. 



(1) 



At time t, we randomly choose a previous time 1 < t' < 
ft with equal a priori probabilities. The walker then 
chooses the current step direction Vt based on the value 
of Vf, using the following rule. With probability p the 
walker repeats the action taken at time t', and with prob- 
ability 1 — p the walker takes a step in the opposite di- 
rection —Vf- Values p > 1/2 (p < 1/2) generate pos- 
itive (negative) feedback. For p sufficiently larger than 
p = 1/2, the behavior becomes persistent (i.e., H > 1/2). 
But the finding of persistence for p < 1/2 and small / 
overturned commonly held beliefs concerning repetitive 
behavior and memory loss Q. Very recently, Kenkre Q 
has found an exact solution to this problem for the be- 
havior of the first moment, for all /, and generalized it in 
important ways, with excellent agreement with the nu- 
merical results over 6 orders of magnitude in time. Here, 
we investigate how persistence depends quantitatively on 
memory loss and how to characterize the important un- 
derlying symmetry properties. 

Fig. 1(a) shows values of H(f,p), estimated via simu- 
lations, as a function of the feedback parameter p and the 
memory fraction /. We choose an order parameter 2H—1 
that has positive values only in the persistent regime. 
Misleadingly, only two phases may at first seem appar- 
ent, namely, persistent and nonpersistent. However, the 
persistent regime comprises two different phases with dis- 
tinct symmetry properties. For p > 1/2 we find classical 
persistence satisfying continuous scale invariance sym- 
metry. In contrast, we find discrete scale invariance sym- 
metry for p < 1/2: scale invariance holds only for 
discrete values of the "zoom" or magnification = A fe 
(k = 1, 2, 3 . . .). The nonpersistent regime, similarly, al- 
lows for monotonically growing as well as log-periodic 



2 



solutions. The simulation results agree well with our an- 
alytical results discussed further below (see Fig. QIb)). 
This spontaneous symmetry breaking indicates a distinct 
phase (and eliminates the possibility of an infinite-order 
phase transition) [lOfl. Mathematically, discrete scale in- 



variance symmetry [J] involves complex fractal 11 . 
dimensions: for z = a + bi, the real part of t z equals 
t a cos [6 log i], indicating log-periodicity [3]. Fig. 1(c) 
parametrically plots two independent persistent walks for 
p < 1/2 (p = 0.1 and / = 0.1). The walks appear not 
only statistically but also geometrically self-similar. In- 
deed, we find a pattern reminiscent of a logarithmic spi- 
ral. We estimate a value of the critical exponent = 1 
from the double log plot (Fig. 1(d)) of the order param- 
eter versus \p — p c \. 

We next approach the problem analytically. We choose 
ty = ±1 and use the previously discussed recurrence 
relation, Eq. [T] Let Ufit) e n&(i) denote the number 
of steps taken in the forward and backward directions 
respectively, at time t (inclusive). The total number of 
steps is thus n/(£) + rib(t) = t. For full memory, the 
probability to take a step in the forward direction at time 
t+1, fori > 1, is P+{t) = ^W^(l-p). Similarly, 

p- s (t) = ^-p + - p). So, the effective value 

expected at time t + 1 is vf +1 = P^(t) — P^g(t)- Since 
Xt — n>f(t) — nb(t) +xq, we obtain vf +1 = a Xt ~ x ° , where 
a = 2p — 1. We can interpret this result as a series of 
experiments or walks at time t having the same number 
of steps forwards and backwards, giving the value vf +1 . 

Now we introduce memory loss. Let the memory range 
be L = L(t) = int(ft), where int(x) denotes the integer 
part of x, for < / < 1, starting at t = 0. In analogy to 
the results above, we obtain vf +1 = a XL J^ Xo . 

Now we study the nth moments of x™. Taking Eq. Q]to 

" / n\ 

power n we get x% +l = (x t +v t+ i) n = J2 I . J v\ +1 x™ \ 

For all even exponents i we know v\ +l = 1 and for odd 
exponents vl +1 = v t +\. Using the expression for vf +1 , 
with xq = 0, we obtain 
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where A = ^ and s(n) 



n-A-l 



We have A = 1 



2 " " "VJ ~ 2 

for even n and A = for odd n. If s(n) < 1, then the 
sum vanishes. 

In the asymptotic limit, we arrive at the following dif- 
ferential equation for the moments, starting from Eq. [21 



~dt { 
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FIG. 1: (a) Hurst exponent H(p,f) estimated from simu- 
lations of non-Markovian walks that remember only a frac- 
tion / of their distant past, with feedback parameter p. 
Color scheme: H < 0.6 (white), 0.6 < H < 0.7 (green), 
0.7 < H < 0.9 (red), H > 0.9 (black). Full circles show 
the edge where H = 1/2. Persistence {H > 1/2) arises for 
positive (p > 1/2) as well as negative (p < 1/2) feedback, 
(b) complete phase diagram showing the 4 phases, plotted 
according to the exact solutions given by Eqs. 1111 1131 The 
dashed line fo(p) delineates the threshold for log-periodicity 
and cleaves the nonpersistent regime into two (Eq. |HJ. (c) 
parametric plot of positions x\(t) and X2(t) for two realiza- 
tions of log-periodic walks. Inset shows zoom of the center 
(boxed area) of the logarithmic-like spiral. Note the geomet- 
ric self-similarity, (d) double log plot of 2H — 1 versus p c — p, 
showing critical behavior with critical exponent (3 = 1. 
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For the first moment (n = 1), we obtain an equation 
identical to the one obtained by Kenkre Q: 



d a 

di {xt) = Tt { ft) - 

Considering an expansion of 



the 



(4) 
form 
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FIG. 2: The order parameter 2H — 1 versus p and /, estimated 
using Eqs. [1(31 and [HI 



{x t ) ~ J2i Ait 6 * sin(B l hx(t) + Cj), we obtain 
tem of transcendental equations for B and 6: 



6 
B 



" 1 cos(Bln/) 
^sin^ In/). 



sys- 

(5) 
(6) 



For a > 0, we have B = 0. The system given by 
Eqs. [5] and [6] reduces to Eq. [7] below, ruling out a log- 
periodic solution. For 5 > 1/2, we obtain superdiffusion 
(see below). 

For a < 0, there exists a threshold denned by a contin- 
uous set of values (p, /), with oscillating solutions. Con- 
sider first the case B = 0, without oscillations. Then, 
Eq. [5] becomes 



af 



s-i 



(7) 



which only has solutions for / > /o (p) . Numerically solv- 
ing this equation, we find the threshold /o(0) = 0.7569 
for p — 0, in perfect agreement with the expression 



(l-2p) ln(l// ) = f /e 



(8) 



obtained by Kenkre (personal communication, 13 July 
2007) for the onset of log-periodicity |3j, shown as a 
dashed line in FigQJb). Indeed, we have found an al- 
ternative proof of Eq. [5] starting from Eq. [71 using the 
Lambert W function. For S > 1/2, both superdiffusion 
and log-periodicity appear (see below). 

We next study the second moment. If n = 2 then 
A = 1 and s(n) = 0. Thus, Eq. [3] leads to 



dt 



<*?> 



2a 



Xf t X t ). 



(9) 



Using the fact that \(x)\ < ((x 2 )) 1 / 2 , we can prove 
that there exists a function A(t) such that (xftXt) = 
A(t)((x%)(xl))V*,with-l<A(t)<l. 

For a > (p > 1/2), no oscillations appear and 
A(t — > oo) = 1. We can thus show that the following 
transcendental relationship holds for the Hurst exponent, 
asymptotically : 



H-l 



H = af 

This result is identical to Eq. [7J with 5 
quently, the curve 



(10) 



H. Conse- 



fo = 16 [Po - 



1 



(11) 



separates the diffusive and anomalous regions for p > 1/2 
in the (p, /) plane (Fig. O^b)). The case / = 1 leads to 
p c = 3/4, in agreement with ref. [l[. 

For a < 0, we try the expansion 
( x t ) ~ Si a it 2Hi sin 2 (6i lni + c,). We assume that 
the dominant terms of ( Xt ) and ( x\ ) have the same 
"period" and phase difference, so that b = B and c = C. 
In the log-periodic region (i.e. 6 ^ 0), we also assume 
that H > 6, so that the solution follows from Eqs. [5]e[Sl 
Indeed we conjecture that for walks lacking subdiffusion, 
6 > 1/2 implies H = S always. Subject to natural 
restrictions, the Hurst exponent must thus satisfy 



H 



2H-2 



H 2 



tan lnC/Va 2 / 2 ^ 2 - H 2 



(12) 



for 5 > 1/2 and H — 1/2 otherwise. The separation 
line of the diffusive and anomalous phases corresponds 
toiJ = <5=l/2, so the critical line (Fig. [TJb)) satisfies 




a 2 1 



(13) 



We obtain the critical value of / c (0) = 0.3284 for the on- 
set of log-periodic superdiffusion, which occurs at p = 0. 
Note that p c — 1/2, f c — represents a multicritical 
point. Fig. [lib) shows the complete phase diagram, con- 
sisting of 4 phases. Fig. [2] shows a better view of how the 
order parameter 2H — 1 depends on p and /. 



Substituting our choice of 
definition of A(t), we obtain 



) into Eq. [9j with the 



H sin(b In t + c) 2 + b sin(6 In t + c) cos(6 In t + c) 
af 11 ' 1 |sin(6 In t + c) \ \sin{b In ft + c) 



A(t) 

u. i " " ibimoui i i- c) i launom /i-|-cj| 

(14) 

Semi-empirically, we find that A(t) ks ±1. Specifically, 
we propose 



A(t) 



sin(6 In t + c) cos(£> In /t + c) 
| sin(& In i + c) sin(6 In /i + c) | 



(15) 
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Note that A(t) = 1 for / = 1. Applying the same rea- 
soning for the case H —1/2, we obtain 

Alt) = / 1/2 [2aa|sin(61nt + c)| |sin(&ln/t + c)| ] _1 
x [a sin(61nt + c) 2 
+2a&sm(&lni + c) cos (feint + c) - l] . (16) 

Thus, exactly on the critical line, we get a marginally 
superdiffusive solution: (x 2 {t)) — at In t sin 2 (b In t + c). 
simulations suggest that higher order terms in the expan- 
sion become important near the critical line. 

We next focus on finding a suitable Fokker- Planck 
Equation (FPE). Let Y denote the position of a walker. 
Consider the conditional probability at position Y at 
time t + 1, given position xq at time t = 0: 

P(Y,t +l\x ,0) = P(Y + l,t\x ,0)P* e (t,Y + 1) 

+ P(Y-l,t\x ,0)P^(t,Y-l). (17) 

Using the definitions of rif(t) and rib(t), and the proba- 
bilities to go forwards or backwards, we obtain 



1 ± 



a(Y-x ) 



(18) 



Substituting Eq. QS in Eq. [T7| we get 



P (Y,t + l|i ,0) 



1 - 



a(Y-x + l) 



a(Y -x Q -l) 



t 



P(Y + l,t\x ,0) 
P(Y-l,t\x ,0). (19) 



From this last equation, in the asymptotic limit for t and 
Y, we get a FPE for x = Y which is identical to the 
one in ref. [l[, as expected. Similarly, we can show for 
memory loss that 



dP(x,t) _l d 2 

Ot ~ 2dx 2 {,) 



fllxftPM], (20) 



where Xf t denotes the position at time ft that leads to 
position x in the future, at time t. 

In summary, we have uncovered the essential features 
of the phase diagram for this problem, based on numer- 
ical as well as analytical results. We expect the phase 
diagram and other findings reported here to contribute 
towards a better quantitative description of persistence 
in diverse economic [l3| , sociological [3 , ecological [l5| , 
biological 0, UH and physiological [^complex systems 



where recent memory loss may play a role Q. Specif- 
ically, prime candidates for further study include sys- 
tems with long-range memory that show evidence of log- 
periodicity and discrete scale invariance symmetry [4|. 
The most important result, for systems with negative 
feedback, concerns the existence of a critical threshold of 
memory loss for the onset of superdiffusion. Finally, we 
hope that the insights provided by this study dispel lin- 
gering doubts that negative feedback coupled with recent 
memory loss can in fact cause persistence. 
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